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Asymptotic freedom of gluons in QCD is obtained in the leading terms of their renormalized 
Hamiltonian in the Fock space, instead of considering virtual Green’s functions or scattering am¬ 
plitudes. Namely, we calculate the three-gluon interaction term in the front-form Hamiltonian for 
effective gluons in the Minkowski space-time using the renormalization group procedure for effective 
particles (RGPEP), with a new generator. The resulting three-gluon vertex is a function of the 
scale parameter, s, that has an interpretation of the size of effective gluons. The corresponding 
Hamiltonian running coupling constant, g\, depending on the associated momentum scale A = 1/s, 
is calculated in the series expansion in powers of go = gxg up to the terms of third order, assum¬ 
ing some small value for go at some large Aq. The result exhibits the same finite sensitivity to 
small-a; regularization as the one obtained in an earlier RGPEP calculation, but the new calculation 
is simpler than the earlier one because of a simpler generator. This result establishes a degree of 
universality for pure-gauge QCD in the RGPEP. 

PACS numbers: ll.lO.Gh, 11.10.Hi, 12.38.-t, 12.38.Aw 


I. INTRODUCTION 

This article describes a calculation of asymptotic freedom m in the leading ultraviolet terms of front form [5] 
(FF) Hamiltonians for gluons, using the renormalization group procedure for effective particles (RGPEP, see Sec|lT| 
developed in recent years as an element of the program of constructing non-perturbative QGD outlined in Ref. [1]. 
Besides the asymptotically free ultraviolet behavior, the calculation also conhrms finite dependence of the effective 
Hamiltonian three-gluon coupling constant on the regularization of small-a: singularities in the bare theory. This 
dependence is of interest since the small-a: behavior may in general be thought related to the vacuum state in the 
instant form (IF) of dynamics and in that form the vacuum is believed to be responsible for symmetry breaking 
and confinement. But the third-order RGPEP calculation reported here would have to be extended to higher orders 
to verify if it can shed any light on the relevant mechanisms. 

In the earlier RGPEP calculation [7] , a generator is used that is suitable for perturbative calculations but difficult to 
use beyond the perturbative regime. In this article, we use a generator that is much easier to use beyond perturbative 
expansion. The difference between the generators is further explained in Sec. |llj At the same time, our perturbative 
calculation demonstrates that the new RGPEP generator passes the test of producing asymptotic freedom, which 
any method aiming at solving QGD must pass. In particular, passing this test is a precondition for tackling non- 
perturbative issues, such as the ones that emerge when one allows effective gluons to have masses [S]. 

An additional point is thus made that two different versions of the RGPEP, defined using two different generators, 
yield the same behavior of the running coupling constant in three-gluon Hamiltonian interaction term when the calcu¬ 
lation is carried out in the third-order expansion in powers of the coupling constant. This means that a considerable 
change in the the RGPEP generator does not influence the finite behavior of the three-gluon term. Therefore, we 
suggest that there exists certain degree of universality in the behavior of EE Hamiltonians in the RGPEP: the leading 
terms in the Hamiltonian beta function are universal and they are universally obtained using different versions of the 
RGPEP. 
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The paper is organized as follows. The RGPEP is generally but briefly described in Sec. The bare Hamiltonian 
for gluons is introduced in Sec. |III[ which includes the derivation of its FF density from the Lagrangian density, solving 
constraints, quantization, and regularization. Sec. |IV| explains our third-order calculation of the effective Hamiltonian 
three-gluon interaction term as a function of the RGPEP scale. This scale is denoted hy t = where A = 1/s is the 
invariant mass width of the form factors in effective Hamiltonian vertices that solve the RGPEP equation, and s is the 
parameter that has the interpretation of size of effective gluons. The calculation includes identification of second-order 
mass counterterms in Sec. |IV A[ third-order counterterm for the three-gluon term in Sec. |IVB[ and the definition and 
result one obtains for the Hamiltonian running coupling constant in Sec. |IV C| Comparison of the previous and current 
calculation is summarized in Sec. |Vj and Sec. |VI| concludes the paper. Several appendices provide details necessary 
for completeness of the paper, including details of the RGPEP in App.fA] details of the initial Hamiltonian in App.[B| 
integration of the RGPEP equation order-by-order in App. calculation of the three-gluon vertex counterterm in 
App.|Dj and method of evaluating the third-order contributions to the running coupling constant gx with a formula 
explaining the infrared stability guaranteed by the design of the RGPEP in App. 


II. THE METHOD OF CALCULATION 

In distinction from calculations of Euclidean Green’s functions w and from early calculations using infinite 
momentum and light-front techniques [HHH], or other approaches, e.g. [iMm, some discussing three-gluon cou¬ 
pling [IHl HH] , asymptotic freedom of gluons is derived here as a feature of the Minkowski space-time FF Hamiltonian 
that acts in the Fock space of virtual gluons obtained as a result of an explicit operator renormalization group trans¬ 
formation. The transformation does not involve any wave function renormalization constant and the counterterms 
are calculated without assuming multiplicative renormalizability. More specifically, we use the RGPEP mentioned 
in Sec. |l] and we calculate the coefhcient in front of the Hamiltonian operator interaction term that annihilates one 
effective gluon and creates two, or annihilates two and creates one. This coefficient is called the coupling constant. It 
is denoted by gt, since it depends on the RGPEP scale parameter t that corresponds to the size s of effective gluons, 
t = s^. In the calculated effective Hamiltonian operator, the gluon interaction vertices are softened by the form factors 
of width A = 1/s in momentum variables. Asymptotic freedom is exhibited in the behavior of the coupling constant 
gt as t approaches zero, or A tends to infinity. 

The size parameter for gluons is introduced by solving the RGPEP differential equation, 

±nt = [Gunt] , ( 1 ) 

where Ht denotes the Hamiltonian of interest and Gt plays the role of a generator of the required transformation, 
see Appendix The transformation changes the bare creation and annihilation operators for point-like gluons of 
canonical QCD, which are denoted by oq in reference to t = 0, to the operators for effective gluons of finite size 
s = that are denoted by at, 


at 


lit ao U\ 


( 2 ) 


where 


Ut 



( 3 ) 


and T denotes ordering in r. The initial condition for solving Eq. Q is provided at t = 0 by the canonical FF 
Hamiltonian with modifications implied by its divergent nature. These modifications include the regularization factors 
in interaction vertices and the counterterms whose structure is found using solutions to Eq. 0. 

Our choice of the generator has the form of a commutator, similar to but different from Wegner’s [5D], 


j 


( 4 ) 


where operators %/ and %pt are defined using parts of the Hamiltonian Ut ■ The operator Hf is the free part of Ht 
and T-Lpt is defined in terms of the interaction terms as explained in Appendix]^ The definition of TLpt secures that 
Ut is invariant with respect to the seven-parameter Poincare subgroup that forms the kinematical symmetry group 
of the FF of Hamiltonian dynamics. The generator Gt in Eq. Q is simpler and more suitable for non-perturbative 
calculations than the one used in Ref. [7] (see below). 

Once the generator is chosen and the counterterms that complete the definition of the initial condition for T-Lt at 
t = 0 are found, the Hamiltonian for effective gluons of size s is uniquely determined, up to the value of the coupling 
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constant go at some arbitrarily chosen value of t = to- Thus, the renormalized FF Hamiltonian for QCD could in 
principle be defined using the RGPEP without any reference to perturbation theory. 

However, our calculation is only carried out using expansion in powers of gt [U] up to the third order. The reason is 
that too little is currently known about the counterterms required for non-perturbative calculations. Even the terms 
of second and third order require study. We show in the next sections that the generator Qt of Eq. Q produces the 
same third-order dependence of gt on t as the one obtained in Ref. [7]. We thus obtain the Minkowski space-time 
Hamiltonian example of the universality of leading perturbative terms in the coupling constant in asymptotically 
free theories. Our perturbative calculation indicates what kind of terms are necessary to counter the ultraviolet 
divergences. They also show how the small-a: singularities appear in addition to the ultraviolet ones, and that they 
cancel out, leaving behind finite effects. 

The current level of knowledge about the EF Hamiltonians for QCD being quite limited, calculations of higher order 
than the third one discussed here are required to gain more information. However, the higher-order calculations require 
the third-order output reported here as an input. Given the RGPEP Eq. Q and its systematic expansion EH, one may 
hope that it will eventually become possible to identify the structure of Ht required for obtaining non-perturbative 
solutions with mathematically estimable precision, perhaps using conceptual analogies between the RGPEP and 
procedures discussed in Refs. EH E3] . 


III. CANONICAL HAMILTONIAN FOR GLUONS 

The initial condition for solving Eq. Q is the canonical FF Hamiltonian for gluons in QCD plus counterterms. In 
this section we describe the canonical Hamiltonian for gluons. The description introduces the notation for details that 
appear throughout the article. 

The canonical Hamiltonian is derived from the standard Lagrangian density 

£ = , (5) 

where F^'' = d^A" — d’'A^ + ig[A>^, A''], A^ = A°-f^t°‘, and = (5“^/2. The associated energy- 

momentum density tensor reads, 

^ . ( 6 ) 

The FF Hamiltonian is obtained by integrating the component T’*’ ~ over the hyperplane defined by the condition 
x'^ = = 0. We work in the gauge A+ = 0, in which the Lagrange equations constrain A~ to 

A- = ^2d^A^-^ig[d"^A-^,A-^], (7) 

so that the only degrees of freedom are the fields A-^. The first term in A~ is independent of the coupling constant 
g. This term is by dehnition included in a new constrained field, which is denoted by the same symbol A in what 
follows, with 


A- = ^2d^A^ , ( 8 ) 

while the second term is explicitly included in the interaction Hamiltonian that is written in terms of fields A -^, using 
A~ defined in Eq. ([^. Employing this convention and freely integrating by parts, one obtains the EE energy of the 
constrained gluon held in the form 



where % = is a sum of four terms, denoted as in Ref. Eli 


= 77^2 -I- 77^3 -f 77^4 + 'H[aAAY ■ 


The terms are [SHIQ] 


(9) 


( 10 ) 


77^2 = 


( 11 ) 


4 


, ( 12 ) 

Um = -\g^A^,ApnA^,APr , (13) 

'H[9aay = [id+A^, A^Y A^, A^f . (14) 

The bare expression for the quantum gluon energy operator is obtained through replacing A^ in “ by an operator 
four-vector A^, which is defined by its Fourier composition on the front corresponding to x'^ = 0, 




E/w 


^ ^ k/T^k(TC^ 


—iki 


, .c M* t 
^ ^ ^ka^kac^ 


ikx 


c+=0 


(15) 


This operator acts in the Fock space spanned by states created by products of the creation operators on the bare 
vacuum state |0). 

In the operator A^, the Fourier-like integral over kinematical momentum variables is carried out with the measure 
[fc] = 0{k^)dk'^cPk^). Thus, the integration matches the one in the Fourier transform only in the transverse 
directions, the integral over k'^ being limited to positive values. 

The polarization four-vectors e have components 


e 


k(7 


i4a = = 2fc-Le^/fc+,£^) . 


(16) 


The symbol a labels the gluon spin polarization and c is a color index. The creation and annihilation operators satisfy 
commutation relations 


^k<7C1 


a 


t 

k'(7'c' 


k+6{k - k') (5“' , 


(17) 


where 6{p) = 16 tt^6{p'^)S{p^)S{p^), and commutators among all as, and among all a^s, vanish. By definition, a/co-c|0) = 
0 for all momenta, spins and colors. 

Normal-ordering of the operator density 'H(A) defines the FF integrand in the Hamiltonian, 

P~ — ^ J dx~(P‘x'^ : T-LiA) : , (18) 

in which all annihilation operators are on the right side of all creation operators. Details of P~ are given in Appendix [ b| 


A. Regularization 

The bare Hamiltonian is regularized by introducing regulating factors, denoted by r, in the interaction terms. These 
factors make the interaction terms vanish [7] if the associated change of any gluon relative transverse momentum were 
to exceed the very large cutoff parameter A. Likewise, the interactions are also made to vanish if any change of any 
longitudinal momentum fraction x of any gluon involved in the interaction were to be smaller than the very si nal l 
cutoff parameter 5. Various regularization factors can be incorporated in the interaction terms in P~ of Eq. (18) 
to realize these conditions. Appendix shows how the regularization factors are introduced in P~ according to the 
following rules. 

In every interaction Hamiltonian term every particle creation and annihilation operator is labeled by its momentum 
quantum numbers p'^ and p'^. Let the total momentum of all quanta annihilated in a term have components P'^ 
and . These are the same as components of the total momentum of quanta created in the term. The relative 
momentum fraction x for the quantum of momentum p is defined as the ratio 

Xp/P = P+/P+ , (19) 

and the relative transverse momentum for the quantum is defined by 

Kp/P = p^ — xP^ . ( 20 ) 

Every creation and annihilation operator in every term in the entire canonical Hamiltonian of any momentum p is 
multiplied by the regulating factor 


rAs{K^,x) = ri^{n-^)rs{x)6{x) . 


(21) 
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We use one transverse regulator factor 

r^iz) = exp(-z/A^) , 

and one of the following three different small-a; regulator factors, 

a) rs{x) = x/{x + 6) , 

b) rs{x) = 9{x - (5) , 

c) rs(x) = x^ 9{x — e) . 


( 22 ) 


(23) 

(24) 

(25) 


Dependence on the transverse regulator factors will be removed using the RGPEP. Effects of the small-x regularization 
will be described by comparing results obtained using different regulator factors in Eqs. (23) to (25). 

Regularization factors in canonical QCD terms that are quartic in gluon field operators are additionally specified 
by treating every such term as built solely from vertices in which one quantum is changed to two or vice versa. This 
regularization choice also applies to the seagull terms that result from constraints on A~, as if the constrained field 
component corresponded to an exchange of a quantum with a corresponding momentum. Details are available in 
Appendix [ b| In an abbreviated notation, the regularization uses symbols rp^p = rA 5 {P,p), where 


XAsiP ^p/pP iXp^p^VAS^P P (1 Xpjp^P ,1 ^p/p] 


(26) 


B. Counterterms 


The initial condition for solving the RGPEP Eq. ([^ is provided by the regulated canonical Hamiltonian plus 
counterterms. The ultra-violet divergent parts of these counterterms, depending on the regularization parameter 
A, are found in a process of calculating Hamiltonians with finite parameter t and eliminating their dependence 
on A by adjusting the initial condition. More precisely, one adjusts the counterterms so that the coefficients of 
products of creation and annihilation operators in an effective theory for gluons of finite size s become independent 
of the regularization parameter A when the regularization in dynamics of gluons of size zero is being removed. The 
remaining unknown finite parts must be adjusted to respect symmetries of the theory and to match its predictions 
with experiments. In the case of pure glue theory, the only unknown parameter would be Agcp, which could be 
adjusted so that, for example, the theory yields the desired value of mass for some glueball, if the mass gap is found 
to exist. 

The question arises is how the small-x regularization effects can be removed. The required counterterms are relevant 
to our understanding of the theory ground state and mechanism of confinement [B] . We shall show in the next sections 
that the small-a: divergences cancel out in the third-order RGPEP coupling constant in the effective Hamiltonians. 
However, the third-order effective three-gluon interaction terms exhibit a finite small-a; regularization dependence 
which is not yet fully undertstood. 

One hopes that the finite small-a; regularization dependence, which we illustrate using different small-a; regulator 
factors listed in Eqs. (23) to (25), may cancel out in the mass eigenvalues for glueballs and their scattering amplitudes. 


Verification of such cancellation is foreseen to be difficult because it involves solving bound-state eigenvalue problems 
for gluons. It may turn out that even for the calculation of lightest glueballs one needs to introduce the counterterms 
that also secure confinement of color [B]. Before this issue is resolved, in the practice of approximate calculations of 
observables, one can seek finite parts of the counterterms using as a guiding rule the minimization of dependence on 
the renormalization scale [14] . which in the RGPEP means minimization of dependence on the gluon size parameter 


We show in the next sections that there exists a small-x regularization that yields the third-order effective coupling 
constant which depends on the size of effective gluons in the same way as the running coupling constant calculated 
using Feynman diagrams for off-shell Green’s functions depends on the virtuality of external gluon lines. Moreover, 
we demonstrate below that two different RGPEP generators lead to the same third-order results for the effective 
Hamiltonian coupling constant including finite effects of the small-a; regularization. These results suggest that the 
calculated finite sensitivity to the small-x regularization is not accidental and, being established here, should be 
further studied as a potentially universal feature of a whole class of FF Hamiltonians for effective gluons. 


IV. CALCULATION OF THE THREE-GLUON TERM 

We solve the RGPEP Eq. 0 pertubatively, expanding Hf in powers of the coupling constant g up to third order, 
Pl^t = + Hii g2 i + Pl2l,g,t + + H‘il,g'^,t + Hl3,g^,t + ^22,g2_t + ^21,g3,t + Pl^l2,g^,t ■ (27) 
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The first subscript lists the numbers of creation and annihilation operators contained in a term, correspondingly. The 
second subscript indicates the order in powers of g, and the third subscript indicates dependence on the parameter t. 
For building intuition, we introduce symbols: for mass terms, which have the first subscript 11; Y for three-gluon 

interaction terms, which have the first subscripts 12 or 21; and X for four-gluon interaction terms, which have the 
first subscript 22. Consequently, the powers of g are explicitly accounted for using the following notation: 


Hn,o ^ E , (28) 

—>■ g "^, (29) 

H2i,g + Hi2^g —>■ gY2i + gYi2 , (30) 

^ 22,32 —>■ g^X 22 , (31) 

Hzi,g'^ + Hi^ g2 —>• -I-g^.=,i3 , (32) 

-f^21,g3 -|- Hi2^g3 —>■ g^Yfi21 + g^Yfil2 ■ (33) 


To be faithful to the difference in notation between Ht and Ht, associated with changing bare to effective gluon 
operators, see Appendix]^ we could also introduce symbols like E instead of E, yt instead of Yt, etc. However, it is 
simpler to remember the difference, and write 

Ht = E + g'^gt + gY2it + gYi2t + g^X22t + g'^'^^.it + -f g^Yh2it + g^Yhi2t ■ (34) 

The initial condition Hamiltonian at t = 0, is written using symbols with the subscript 0 in place of t, 

Hq = E + g^gQ + gY2io + gYi2o + g^X 220 + + ff^Si3o -I- 3 ^ 1^210 + g^Yhi2o ■ (35) 


The counterterms that need to be found are included in Hq. 

Besides using the parameter t and its initial value t = 0, we also use the parameter A = whose initial 

value is 00 . The parameter A has the interpretation of momentum-space width of the form factors that appear in 
solutions for Ht- 

The RGPEP Eq. Q for the Hamiltonian (34) reads 


g'^dtgt + gY 2 it + gYl2t + 9^^221 
2 - 


■ 9^%it 




+ 9 Yh2it + 9 Yhi2t 

= [ [E, gYiiPt + gYi 2 Pt + 9 ^Y 22 Pt + 9^^3iPt + 9^^i3Pt + g^Yh 2 iPt + 9 ^Yhi 2 Pt\ , 
X E + g^fit + 5^214 -I- gYi 2 t + 9 ^Xt -\- g^’B.^it + g^'^iQt + g^Yh 2 it + 9^Yhi2t\ ■ 


(36) 


We solve Eq. (36) order-by-order in series of powers of g, which eventually is translated into a series expansion in 


powers of gt ■ Mass squared terms are of the second order and the gluon vertex is made of terms of the first and third 
order. After removing powers of g from the equations. 




X'22t 


' 21t 
“31t 


12t 

“13t 


Yin + Yi2t = [[E, Y2iPt + Fi2Pt], E] , 

[[E, X22Pt + '^31Pt + Si3Pt] , E] 

[[A, l21Pt + Fl2Pt] , ^21* + Fi2t] , 

[[E, Yh2ipt + Yhi2Pt ], E] 

+ [[^:A 22 Pt-|- ^31Pt + ‘^13Pt] : ^ 21 * + Yi2t] 

+ [[E, Y2ipt + Tl2Pt] , gt + X22t + Ssit -|- 5i3t] 


YlnAt + YL^t — 


(37) 

(38) 

(39) 


The running of the Hamiltonian coupling constant gt is encoded in the operator Eq. (39). Solving Eq. (39) requires 


knowledge of solutions to the operator Eqs. (37) and (38). The gluon mass squared counterterm is obtained from the 


second-order equations, and the three-gluon vertex counterterm from the third-order equations. 
The solution for 


Yt — g{Yi2t + ^ 21 *) + g^{Yh2it + Yhi2t) 


(40) 


is written in terms of the bare creation and annihilation operators for canonical gluons and powers of the bare coupling 
constant. The last step in the RGPEP is the replacement of the bare gluon operators by the effective ones at scale t 
and expressing g in terms of gt- 

Details of solving Eqs. (37) to (39) are described in Appendix [C| Here we list the results. The first-order terms are 


the same as in Ref. [7] , see Appendix G 1 
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A. Mass squared term and its counterterm 


The second-order mass squared term for effective gluons that solves Eq. (38), has the form 


A? = ^ , 


(41) 


where the only element that depends on the scale t is the parameter The result for it reads 

2 / dxrs^{x) 


2 2,9 

Mi — Mi + 


( 4^)2 


p 1 pOO 

/ dxrs^{x)y\Yi 2 k\^f dz exp (-2te^) , 
</ 0 22 " ^ 


(42) 


where 


y\Yi 2 k\Vn^ = N,[l + l/x^ + l/{l-xf] = P{x)l[2x{l-x)] 


(43) 
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and P{x) is the Altarelli-Parisi gluon splitting function Pgg(x) [21] ■ Nc = 3 denotes the number of colors. The 
effective mass squared term is sensitive to the small-x regularization. The counterterm that canceled dependence on 
the ultraviolet cutoff A —> oo contains the mass-squared factor of the form 


2 2 , 
Mo — Mi+ 
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[Att)'- 


pi pOO 

/ dxrs^{x)P{x) / di: exp [—42 x(1 — x)/A^] . 
do do 


(44) 


Comparison with Ref. [7] shows that the gluon mass-squared term obtained using the RGPEP generator of Eq. (A4) 
does not differ from the one obtained using the generator of Eq. (A5). 


The ultra-violet finite part of the mass counterterm, depends on the small-x regularization parameter 5 in 
the initial Hamiltonian. Therefore, the simplest way of choosing the ultra-violet hnite part of the mass-squared 
counterterm is to set the mass squared for effective gluons at some value of t to a desired function of <5. Such function 
of 5 can be fixed by demanding that the effective Hamiltonian eigenvalue for lightest states with color quantum 
numbers of a single gluon contains a specified mass-squared term that depends on the parameter d in a specihc way. 

The right dependence for defining a complete theory of gluons is currently unknown but it is also currently not 
excluded that one can attempt to describe conhnement of gluons by demanding that the gluon mass eigenvalue 
diverges in the limit d —> 0. Verification of this option requires studies far beyond the scope of this article. Namely, 
one needs to study terms of higher order and consider the eigenvalue problem in higher order than third, before 
one will know if the perturbative expansion of the RGPEP can lead to establishment of a general structure of 
the Hamiltonian that may solve Eq. 0 beyond perturbation theory. Here, we shall find that, once the ultraviolet 
divergences are removed, the small-x divergences do not appear in the third-order asymptotically free coupling constant 

in renormalized Hamiltonians for effective gluons. _ 

The second-order effective gluon mass term, denoted by Wj, is obtained from in Eq. (41) by applying the 


transformation lAt and thus replacing the creation and annihilation operators for bare gluons by the ones for effective 
gluons of size s. Namely, 


ml 


= Ut A? u, 


This transformation amounts to the replacement in Eq. (41) of a\^^akcrc 
spond to thin and the latter to thick lines in Fig. 


by 


t k<7C 


(45) 

atko-c- The former operators corre- 


B. Third-order three-gluon term and its counterterm 


We focus our attention on the term 7721,in Eq. (27), knowing that i7i2,g3,t is its Hermitian conjugate. The term 
has the structure 


where 


H2l,g^,t — Ut 7i21 Idl , 


7t21 — ft 7t21(n) • 

n 


(46) 

(47) 










The factor ft in front of an operator means that the vertex functions in the operator are multiplied by the form factor 


defined in Eq. (A 6 ). 


The subscript n in the sum ranges over ten values, denoted by alphabet letters from a to j. Each of the summed 
terms results from some specific operator product in Eq. (39). In each of these terms there appears a vertex function, 
denoted by 7 („), in the otherwise universal pattern of the formula 


7t21(n) 


123 • 


[123] S{ki + k2 — ks) 
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1 


2 


7(n) U 1 U 2 U 3 


(48) 


The vertex functions 7 (n)i with subscript n ranging from a to j, are given in Appendix |C 31 with their operator 
origin in Eq. (39) being illustrated by diagrams in Fig. there. The thick external lines in h'lg. correspond to the 
creation and annihilation operators that appear in the three-gluon interaction term for effective gluons of size s. The 
thin internal lines correspond to the commutators that result from moving annihilation operators to the right of all 
creation operators for gluons of size zero, in terms of which the RGPEP Eq. 0 is solved. Lines with a transverse 
dash indicate instantaneous interactions in Hq that result from the constraint of Eq. ([^. 

The vertex functions diverge when the ultra-violet cutoff parameter A is being sent to infinity. The divergences 
result from integration over the transverse relative momentum of virtual quanta whose creation and annihilation 
operators were contracted in the products that appear on the right-hand side of Eq. (39), cf. Eq. (C 8 ). The required 
counterterm, calculated in Appendix has the form 


Yh2io = I [123] Hki +k2- fcs) 

123 


IGtt^ 


7o oja^aa 


(49) 


with the vertex function 70 derived in Eq. (D12), 


70 = -Yi23^ In — {fVc [11 + ^(a^l)]} + 7 finite • 

o IJj 


(50) 


For calculation of the Hamiltonian running coupling constant, the finite part of the counterterm, denoted in Eq. (50) 


by 7 finite 5 will not need to be specified when the subtraction of the diverging part is introduced as described in the 
next section. 


C. Running coupling constant 


Our Hamiltonian running coupling constant gt is extracted from the three-gluon terms in Hf. These terms create 
two gluons and annihilate one, or vice versa. Both types yield the same result for gt- The three-gluon term is a 
sum of terms denoted by (a) to (j) in the previous section. The vertex function of the entire sum depends on the 
gluon colors, polarizations and momenta. In the terms that vary with the scale parameter t, the dependen ce o n color 
and polarization in the limit Kj ^2 0 takes the form of a combination yi 23 shown in Appendix [b| in Eq. (B3). This 
combination is multiplied by a function g(t, cci), where Xi = 1 — X 2 refers to the -f-momentum fraction carried by one 
gluon of the total momentum of two gluons that are created or annihilated by the three-gluon term. The coupling 
constant gt is defined as the value of g{t,Xi) at some value of Xi = Xq, 


gt = g{t,xo) ■ (51) 

Note that the Hamiltonian Tit that appears in Eq. ([^ is calculated using the bare creation and annihilation operators 
and the effective Hamiltonian Ht is obtained from Ht by inserting the effective creation and annihilation operators in 
place of the bare ones. The vertex function is not changed. We calculate gt using Hf 
The counterterm contribution to the vertex function can be written as 


70 7to A 7finite 


(52) 


The diverging part of the counterterm is thus specified using the negative of 74 at an arbitrary finite value of to- 
Therefore, yonite may differ from 7 finite in Eq- (50) by terms that do not depend on t. Such terms will not contribute 
to the dependence of gt on t and will not be further discussed in this article. Note that since the diverging part of the 
counterterm may be a function of xi, as displayed in Eq. (50), one also has to consider the counterterm finite part 
that may be a function of xi [ 6 j. 
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After inclusion of the counterterm defined in Eq. (52), our result for the three-gluon interaction term in "Ht has the 
form (the symbol a stands for spin variables) 


Yt = gYrtYg^YztY ... , 

Yt = ^ Ai 23] 5(1-h 2-3) fi 2 t b)(a;i,K^ 2 >o-) al4«3 , 

123 

Yu = ^ /[123] 4i + 2-3) /i2t Yitixi, Ki 2 ^ 44a3 , 

123 

Yst = ^ Ai23] 4i + 2-3) fi2t Y3tixi,Ki2,(T) a\a\a3 . 

1 OQ 


(53) 

(54) 

(55) 

(56) 


The symbols Y denote vertex functions without the form factor ft. Assuming a counterterm that involves a subtraction 
at some t = to as described above, one obtains the third-order vertex factor Y of the structure 


Y3t{xi,Ki2^<^) = T3t{xi,K32,Cr)-f3to{xi,K^2^a)+f3finite{xi,Ki2^>^)^ 
where the symbol T denotes the sum of third-order terms from (a) to (i) calculated in Appendix 

i 

T = • 


C3 


(57) 


(58) 


The diverging part of the counterterm denoted as 7 (j) in Eq. (D14| in App. C3 is included here through the 
subtraction at t = to, and the associated change in the finite part is not needed iii the discussion that follows. 
Combined, all the three-gluon terms in expansion up to third-order in powers of g, have the form 


Yt{xi,Kj- 2 ,a) = 5Yit(a:i,K4.0')+4 Tstixi, K32,cr) - fstoixi, Ki 2 ,cr) + fsdniteixi, K 32 , cr) 


(59) 


By our definition, the Hamiltonian coupling constant gt i s fou nd as a coefficient in front of the canonical color, spin 
and momentum dependent factor Yi 23 (a;i, k]*‘ 2 , cr) of Eq. (B3), in the limit Kf 2 0, for some value of Xi, denoted 


by xq. At some arbitrary value ot t = to, gto must be set to a specific finite value go that produces agreement with 
experiment when one describes data using the Hamiltonian with t = to. 

We obtain 


lim Yt{Xi,Kf 2 ,Cr) = lim Ct{xi, Kf2)Yl23ixi, Kf 2 ^<^) + yYsHniteiXl, 

KjL-i-O L 

= lim gYi23ixi,Ki2><^) 

-H lim g"^ [c3t{xi,K32) - C3to{xi,K32)\Yi23{xi,K32,(^) 


(60) 


-I- lim g^ T 3 


9—>-0 


3 finite (^1 5 ^12 ) ^ 


) ■ 


(61) 


Removing Ti23(a:i, ^ 12 ;'^) from all terms besides the term T 3 that does not have to have the spin and momentum 
structure of 1123 ( 3 ^ 1 , ^^- 2 , a), in the limit, 


or 


lim ct{xi,K^ 2 ) =5 + 4 lim [c3t(a;i, ^4) - C3to(a;i, k^)] : 


ct(a;i,0-^) = g + g^ [c 3 t{xi, 0 -^) - C 3 ta{xi, 0 ^)] . 


(62) 


(63) 


Evaluation of the limit in Eq. (62) that yields Eq. 


IS 


described in Appendix ^ 


Setting Xi = xo and dropping the argument k ^2 set to zero, our definition of t 


le coupling constant gt leads to 


gt = Ct{xo) 

= g + g^ [cstixo) - C3to{xo)] . 


(64) 

(65) 
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We calculate the coupling constant g that appears in the initial Hamiltonian after inclusion of the counterterm, by 
demanding that at t = to the coupling constant should have the value go, 

9to = 90 ■ (66) 


The value of go is determined by comparison with data using the Hamiltonian corresponding to t = to- Hence, with 
accuracy to terms of order g^ or smaller, Eq. (65) implies 


9t =90+ gl [c3t{xo) - C3to(xo)] 


(67) 


The right-hand side of this result is calculated in Appendix]^ It can also be expressed as a function of momentum 
scale A = 1/s, which facilitates comparison with Refs. [312] and |7]- For this purpose, we denote gt by g\ when we 
set t = A”'^. 


From all terms that contribute to the right-hand side of Eq. (67), listed as (a) to (i) in Appendix]^ only the 
contributions of 1(d) and 7 (g) are different from zero. These terms yield 


9\ 


[11 + M^o)] In^ 


( 68 ) 


where 


h{xo) 

x{xo) 


A(a;o) + x(l - a:o) , 

6 dx — x)-\-l/{x — xo) + l/x] 

J Xq 


^ fs{xo) / dxrs^ix) 



(69) 

(70) 


This result depends in a finite way on our regularization of small-x divergences. 

For the choices of small-a; regularization that are listed in Sec. HI A in Eqs. 
as versions a), b) and c), the limit <5 —>■ 0 yields the Hamiltonian running coupling with the function h{x) given by, 
correspondingly. 


(24) and (25), to which we refer 


a) 

b) 

c) 


h{xo) 

h{xo) 

h{xo) 


, „ „ ^ — Xo — Xn 

12 3-f , , ° , 

(1 - a;o)(l - 2a;o) 

12 In min(xo, 1 — Xo) , 

0 . 


In Xq 


(1 - Xq)'^ - Xq 
a:o(l - 2xo) 


ln(l - Xq) 


(71) 

(72) 

(73) 


IS 


The running of the Hamiltonian coupling constant described by Eq. ( |68[ ) for small-a; regularizations a) and b) 
illustrated by the dashed curves in plots a) and b) of Fig. respectively. Different dashed curves correspond to 
different values of Xq. Only examples with xq between 0.1 and 0.5 are plotted, because the function H^xq) in Eq. (68) 
is symmetric with respect to the change a:o —>■ 1 — a;o. For the regularization c), we have h(xo) = 0 irrespective of the 
value of Xo, and 


9\ = 9o 


In^ 

487r2 Ao 


(74) 


which is illustrated by one and the same continuous line in both plots a) and b) of Fig. Differentiation of Eq. (74) 
with respect to A produces 


where 



Podl , 


Po = 


UN, 

487r^ 


(75) 


(76) 


This result matches the asymptotic freedom result in Refs. [DE], when one identifies A with the momentum scale of 
external gluon lines in Feynman diagrams. Discussion of this result is provided below in Sec. |V| 
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FIG. 1: The FF Hamiltonian third-order RGPEP running coupling constant for effective gluons, g\ of Eq. ( 681 , is drawn using 
different dashed lines for different values of so, as a function of A in GeV, sta rtin g from an arbitrarily chosen value go = 1.1 at 
Ao = 100 GeV. Plots a) and b) correspond to small-s regularizations in Eqs. (231 and ( |24[ ). The thick continuous lines in both 
plots show one and the same result for the regularization in Eq. (251, which exhibits no dependence on xo- 


V. UNIVERSALITY OF THE RGPEP SOLUTION 


As a result of the third-order RGPEP, the three-gluon interaction term in effective FF Hamiltonians for gluons is 
fftAs = -p) fi 2 t alj^al^ats + H.c. , (77) 

123 

where with accuracy to terms order one has 

Yt{xi,K^2^a) = goYi23{Xi,Ki2,(^) + gl T3t(a;i, ^^2: o') “ ^3to o’) • (78) 

For infinitesimal Ki 2 , 


Yt{xi,Ki2,(^) = Vt{xi)Yi2z{xi,Ki2,(^)+glTmnite{Xl,Ki2,<y)+0{Ki2) 


(79) 


where 


Vtixi) = gt + gt [czt{Xl) - Cst{xo) - Cztoixi) + Cstaixo)] 


(80) 


gt is given in Eq. (67) and the coefficients C 3 are described in Sec. IV C| Universality of this result is claimed on the 
basis of comparison with results obtained in Refs. [Big and [7]. 

Comparison with Refs. mm shows that the FF Hamiltonian running coupling constant g\ exhibits, in the RGPEP 
of third order, the same leading dependence on the momentum width of vertex form factor A, as the running coupling 
constant in Refs. mm exhibits as a function of the length A of Euclidean momenta of external gluon lines in the 
three-point effective action. In order to compare these two results, one has to assume that the Euclidean Green’s 
functions correspond, by some continuation procedure from imaginary to real time variable, to a Minkowski space-time 
quantum theory in which a renormalized Hamiltonian has a three-gluon interaction term of a specific dependence on 
the momentum scale parameter A. Our calculation suggests, but does not prove, that the Euclidean scale A corresponds 
to the RGPEP width A. Namely, the FF Hamiltonian matrix element that appears in the virtual transition amplitude 
between one- and two-gluon states in the Fock space of effective gluons, is suggested to correspond to the continuation 
of the three-point Eculidean Green’s function, or effective action, to the Minkowski variables. The suggestion is not 
verifiable by any simple continuation because we do not fully know the analytic structure of either function. However, 
the observed universality of asymptotic freedom in both the perturbative Euclidean Green’s function calculus and 
Minkowskian Hamiltonian quantum mechanical operator calculus points out a direction in which one can seek a 
constructive demonstration that these two ways of defining a theory are equivalent. 

Comparison with Ref. [7], where the Hamiltonian three-gluon vertex is calculated as a function of the momentum 
scale A using the RGPEP with a different generator than the one used here, is facilitated by observing that the size 
of effective gluons, s, is equal to the inverse of A. Using this relation, one sees that the present result is the same 
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as in [7], despite that the generators are different. Specifically, using Qt in Eq. (A5) instead of the one in Eq. (A4) 
does not change the third-order results for g\. Fig. [^illustrates this finding by showing the present results for g\. 
The current calculation explicitly extends the universality of leading perturbative terms in the beta-function to the 
RGPEP calculus for Hamiltonian operators in the effective particle Fock space. 

Thus, the universality we claim is two-fold. One universal aspect is that the third-order RGPEP Hamiltonian 
running coupling constants exhibit the same asymptotic freedom behavior that is known to be universal in the calculus 
based on the renormalized Feynman diagrams. This is of interest from the point of view that the Hamiltonian quantum 
mechanics in the Minkowski space-time and the Feynman diagrams for virtual transition amplitudes can be precisely 
related to each other in a relativistic theory including renormalization, which generally remains to be desired |22112dj . 
Such relation is needed for incorporating non-perturbative features of hadrons in calculations that so far remain 
limited to the usage of qualitative and quantitative input from the parton model. 

The other universal aspect is that the third-order Hamiltonian running coupling constant depends on the size of 
effective gluons, or momentum width of effective Hamiltonian three-gluon vertex, in a way that does not depend on 
the choice of the RGPEP generator. This is of great interest in view of the fact that the presently used generator does 
not depend on the derivative of the Hamiltonian with respect to t, while the previously used generator does. Hence, 
the obtained stability of asymptotically free behavior of effective gluon interactions, with respect to change of the 
RGPEP generator, suggests a viable way around the difficult problem of solving for the derivative of the Hamiltonian 
in terms of the Hamiltonian itself. The generator used here is thus shown to offer a way of seeking non-perturbative 
solutions to the RGPEP equation in a greatly simplified setup in comparison with the original one. 


VI. CONCLUSION 

Knowledge of the third-order RGPEP result for the Hamiltonian of effective gluons is not sufficient for setting up 
any physical eigenvalue problem, such as the eigenvalue problem for a glueball. At least fourth-order terms are needed, 
which describe interactions among two effective gluons including the effect of running of the coupling constant. Such 
calculations are considerably more involved than the third-order calculations described here. 

However, the presently used RGPEP generator, demonstrated here to imply the Hamiltonian running coupling 
constant of the form that is familiar from other formalisms and renormalization schemes, turns out to lead to a 
considerably simpler third-order calculation than the previously used generator did. The consequence of this result is 
that the required fourth-order calculations with the presently used generator are expected to be considerably simpler 
than they could have been with the previously used generator. 

In particular, the same finite effects of small-cc regularization are found using the present generator and the previous 
one. Since the FF Hamiltonian mechanisms of confinement and chiral symmetry breaking are expected to be related 
to the gluon dynamics at small-x, the simpler generator than the one used before is welcome as a tool for studying 
the small-x dynamics in full QGD. 
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Appendix A: Details of the RGPEP 


The effective Hamiltonian is related to the regulated canonical one with counterterms by the condition of no 
dependence on the arbitrary RGPEP scale parameter t, 


This condition implies via Eq. (j^ that 


— 'Ho{ao) . 


ntiao) =Ulnoiao)Ut . 


(Al) 


(A2) 


Differentiation of Eq. (A2) with respect to t yields 


7^t(“o) = [Gtiao),'Htiao)] , 


(A3) 
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where Qt = i® called the RGPEP generator and the related solution for Ut is given in Eq. 

We consider two different generators, one from Ref. [7], 

Gt = , (A4) 

and another one from Ref. [2T], 


Gt = 


(A5) 


The curly bracket in Eq. (A4) indicates that Gt satisfies the equation [Gt,T~l.f] = (1 — which is desi gned 

according to the similarity renormalization group procedure described in Ref. |25| . The form factor ft in Eq. (A4) 
is chosen in the form that also appears in lowest-order solutions obtained using the generator defined in Eq. (A51. 


Namely, in an interaction Hamiltonian term in which R and L refer to the effective particles that enter and emerge 
from the interaction, respectively, the form factor is 


ft = 




(A6) 


where Mp and Aip denote the free invariant masses of the corresponding particles. 

The operator 'H/ in Eq. (A5) is called the free Hamiltonian. It is the part of 77o(ao) that does not depend on the 
coupling constants. 






(A7) 


where i denotes the quantum numbers of gluons and is the free FF energy for the gluon kinematical momentum 
components and pf , 


P^ = 


Pf 


(AS) 


The operator Hpt is defined in terms of "H*, the latter considered an arbitrary series of normal-ordered powers of the 
creation and annihilation operators. 


OO 

■^((oo) = X! X! ct(ii, ...,i„) • • • aoi„ . (A9) 

n =2 ,22 

Namely, Rpt differs from Rt only by multiplication of each and every term in it by a square of a total -I- momentum 
involved in a term. 


n —2 2l,22v)*n 

This multiplication implies that Rt for all values of t possesses seven kinematical symmetries of the FF Hamiltomian 
dynamics: three translations within the front, rotation around z-axis, two transformations generated by -I- 
and — J^, and the boost generated by . The latter is the seventh symmetry generator that does not have a 
counterpart in the commonly used instant form of dynamics, which has only six kinematical symmetries. 

Solutions to the RGPEP equation can be found using expansion in powers of the coupling constant g. Such expansion 
is used in Sec. EYi The last step in the RGPEP is the replacement of bare creation and annihilation operators by 
effective ones. One obtains 

Ht = UtRt{at,)ul . (All) 

This operator can be used for approximate, i.e., neglecting quarks, computations of the states of hadrons made of 
gluons and the transition amplitudes for scattering, decay and production processes that involve such hadrons. The 
arbitrary parameter t can be adjusted in order to reduce the complexity of any such calculation to minimum. In 
practice, it means that the size of gluons s is chosen to match the inverse of the momentum scale that characterizes 
a process of interest. 

Since solutions to the RGPEP equations involve interactions that are smoothed by form factors, the procedure is 
thought to provide a means for the understanding of the connection between quantum field theory and phenomenolog¬ 
ical models. For example, we need to find the mathematical connection between QCD and the properties of hadrons. 
This problem needs solution irrespective of the form of dynamics one uses, e.g. see mUZ]. 


OO 

'Hptiao) = ^ ^ ct(ii, ...,■*„) 


= l / 


^Oii 


^Oi 


(AlO) 








14 


Appendix B: Details of the initial Hamiltonian 


All interaction terms in the operator P of Eq. (18 1 , are regulated as described in Sec. 


Ill A 


The regularization does 


not change the field operator and the associated free part H/ of the Hamiltonian, besides the initial condition that 
^ 0. The latter condition eliminates the terms that contain only annihilation or only creation operators. If 
they were present in P~, it would produce non-normalizable states by acting on the vacuum |0) and all other states 
built using action of creation operators on | 0 ) [ 2 H] . 

The regularized Hamiltonian terms corresponding to densities given in Eqs. ( [IT| ) to ( [T^ , are listed below in the 
same order. 


r f^± 2 


Ha^ = 


123 


- p) rAs{^, 1) 9 ^123 a\ala 3 + g a^aaOi 


where 


^123 = if 


— „■ fCiC2C3 


^1^2 * “^3^ — ^1^3 ’ 


X2/3 


E 2 S 3 ■ £if 


a^l/3 


(HI) 

(B2) 


(B3) 


with £ = £■*■ and k = Symbols and p denote the total momenta of created and annihilated particles, 

respectively. 




FIG. 2: Bare three-gluon vertex. 


Ha^ = 


1234 



-A'l I 234 O 1 02^3^4 + I2340'\a\a3a4 


t 

1234 ^ 4^30201 


(B4) 


1234 = 3[ri+2.ir4.3 • £^£4 - £^£4 ' £;£D ^ fac,c,fac,c. 


+ r3 + 2 . 3 £ 4.i (£^£5 • £^£4 - £^£4 • £^£ l ) /- 3C2 jac , c 4 ] _ 


(B5) 


A'a'I 1234 — £i+2,i£3-I-4,3 (£*£3 • £2£4 ~ £i£4 • £^£ 3 ) 

+ [^3 ,i£2.4 + £1,3^4,2] (£^£2 • £3£4 “ £^£4 ' £^£3) /“^i^3yac2C4 

+ [^3,2^4 + £2,3^4,!] (£^£2 • £3£4 “ £t£3 ’ £^£4) 


(B 6 ) 


H[dAA]2 = Y / [1234] 12340 !0^0^04-b h.c.)-f 123401020304 . (B7) 

1234 


- _ If™ ^ „*„*„*„ (^1 ^ 2 ) ( 0^3 “b 3 ^ 4 ) j.aci C2 £ac^Cd. 

-[aAA]2 1234 — — g ri-l-2,l?^4,3 £i£ 2 ■ £3£4 ■--b X2)^- - 


I - - „*_* (^1 Xz){x2+X4) „aciC3jac2C4 

+ ?^l-|-3,l£4.2 £i£3 • £2£4 3, 7 / 


+ £3-12,3£4.i£3£2 • £i£4 


{Xi +X 3 Y 
{X3 - a^2)(a:i -b X 4 ) 
(a;3 + 0 : 2)2 


j;ac3C2 ^aciC4j 


(B 8 ) 
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X, 


1 . 


(a;i X 2 ){x^ X 4 ) „aciC2 faCiCi 


[aAAYl23A — ^ri+2,l?’3+4,3 £ 1^2 ■ £324 + a; )2 

(li + X^){X2 + Xa) faC2C4 


~ r3,1^2,4 + ?'l,3?’4.2j £i£3 ■ £ 2^4 
- [^ 3 . 2^,4 + ^2.3^4.l] £i£4 ’ £ 2^3 


{X2 - Xa)"^ 

{X 2 + Xz){xi + Xa) 


^aciC2 j^a 

j^acics ja 
jac\c^ ^ac2C3j 


(B9) 


(a:i - xaY 

In all these formulas, the dot • is used merely to visually separate factors comprised of scalar products of transverse 



FIG. 3: Graphical representation of terms (B5l and (B 81 
1 






FIG. 4: Graphical representation of terms (|B 6 [) and 



polarization vectors. 


Appendix C: Integration of the RG equations order by order 


The expansion in a series of powers of the coupling constant g is inserted in Eq. ([Rand s olve d for the hrst four 
terms, i.e., including terms of order 1, g, g"^ and g^. We use notation adopted in Eqs. p7| to (39). 


1. First order terms 


Integration of the terms order g yields 


Y21t + Yi2t = ft [5^210 + Yi2o] , 


where the form factor /* is given in Eq. (A61. The corresponding Hamiltonian term is 

-p) ft rs{xi) 


H, 


( 1 ) 


= E 


123 


9 Yi 23 al^al2at3 + g Yf^^ a\^at 2 ati 


(Cl) 


(C2) 


Note the absence of ultra-violet regularization factors and the presence of small-cc regularization factor fs{xi). The 
reason is that the form factor ft removes sensitivity to transverse momenta much larger than 1/s but, for massless 
gluons, does not regulate small-a; divergences. The creation and annihilation operators correspond to the scale 
parameter t. 
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FIG. 5: First-order term of the three-gluon vertex. 


2. Second order terms 


Solutions for second-order terms are 

At = Ao + / [[^J/i-(^21P0 + ^12Po)] ,/r(^210 + ^12o)]„ , 

Jo 

X22t = f dr/.^ ^ [[S,/T-(y21P0 + ^12Po)] )/T(i^210 + ^12 o)]x22 > 

Jo 

Saif = /fSaiQ-h/f f dp/i-^ [[£’)/rd2iPo])/Ti^2io]H ) 

Jo 

Sl3f = /fSi30-|-/f f dr/^ ^ [[iii,/ t-121Po] , /Ti^21o]Hj3 • 

Jo 

The gluon-mass term consists of the product of two bare vertices, see Fig. 

At = [^ 120 ^ 210 ]^ • 

The subscript /r indicates that one extracts the mass squared term from the product of operators in the bracket 
leads to Eq. (^. 



FIG. 6: Graphical representation of the second-order RGPEP contribution to the Hamiltonian effective gluon mass 
Thin internal lines correspond to intermediate bare gluons and thick external to the effective gluons. 


3. Third-order terms 


The third-order term needed for evaluation of the Hamiltonian running coupling is 

Yh21t = ftYh210 

+ ft [ [[£;,X22Pr] /r(^21o)]y^2i “ /^M/r (^12o) 

Jo Jo 

+ ft dr f-^ [[EJriY21Po)],tll]Y^^^- ft dTf-^[X22r[EJriY21Po)]]Y^^^ 

+ ft f dr f-^ [[E, fr{Yi2P0)]^31r]Y^^^ • 

Jo 


In an abbreviated notation, 


Yh2it — ft ^ It: 


where n ranges from a to i and n = j for the vertex counterterm. One has 


7t2i(n) = X! / + k2- ks) 2 ">'(") ®i4a3 

1 OQ ^ 


123 ' 


(C3) 

(C4) 

(C5) 

(C6) 

(C7) 
. This 

term. 


(C8) 

(C9) 

(CIO) 
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We list below results for the vertex functions 7 („) for all values of n from a to i. The counterterm is described 
afterwards in the next section. 


a. Vertex function 7 (a) 




r dxrst{x) 

Ixi - x){x - Xi) 




rAt{^ 


^i(a) 

k. 


J ijk 


+ 2 ^(a) 


+ (1^2), 


where 


rst{x) = rs{x)rs{l - x)rs{xi/x)rs[{x - xi)/x]rs[{x - xi)/x 2 ]rs[il - x)/x 2 ] , 

= exp [-2(^4^ + ^ 


(Cll) 


(C12) 

(C13) 


and 


Jjk 


*j 


S(a) - ^2 £3 


1 - 


X 


-tti t2 C 3 




x — Xi 

1 


X — Xi 
-X2 


Xi 

1 

1 — X 


2x X XX 2 XX 2 

-1-1-1- 

Xi {1 — x)xi {1 — x)xi {x — Xi)xi 


1 ^*k *j i 

-X2 

I *k j 

X2 

-t- ti t2 C 3 

1 

I 

ti t2 C 3 

I 

7? 

I 


XX2 


I I s:'ik 4 

+6162 1 0 £ 


X — Xi {1 — x)(x — Xi) 

+ d^^e' 


X2 

^ (1 — x)^ 


3 /I +0 ^3 


*3 ^*k i 
■ t-i t2 £3 

k 


-X2 


XX2 


{1 — x)xi {1 — x)(x — Xi)xi 


+£163 

+£ 2^3 




' (1 — x)a 
1 


f XX 2 

X2 Y 

7? 

I 

1 —X Y 


(1 — a:)(a: — a;i) 


- - - <5*%: 


*k 


XX2 


d^^e 


*k 


-X2 


{x - XiY 


+ d^'^e*. 


X2 


(1 — x){x — Xi )2 
XX2 


x{x — Xi) 


_ Xik *3 _ ^^2 _ 

^ (1 — xY{x — Xi) 


X2 


- d^'^F*^ 

^ ' (1 — a;)^a;i ' (x —’ 


(CM) 


B; 


i(a) 




xMfg - . 2 ka2\ 

(x2Wlgg + WIm) 




fief fes/fi 2 — 1 


Mie + M^ + Mi,-Mi2 


feefbd/fi^ — 1 


X2Mlg + xM?i 
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Mlg + Mfg-(M^-M!2)^ 


{2M^-Ml2) 


f fesfie/fi 2 — 1 


12 . 

fcaf / f 12 — 1 


M^-M\2 + M%g + M\^ 2M'^{M‘^ - MI 2 ) 


(C15) 


b. Vertex function 7 (f,) 


7(6) = 2 ^ J VAtin^ ^ £(6) + (1^2) 


(C16) 


where 


£(6) = £f6)^'^ - £i£ 2 • £3K (^1 - S(6) - ^ + ^*^3 • £2« + ^2£3 • £tK 

with S(t,) = (xi + x)(x 2 + 1 — x)/(x — xi)^ and 


B, 


k+ 2M?{M^ - 


2Mf- -M\2 


(CIS) 
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c. Vertex function 7 (c) 


where 


7(c) = 2 


dxrstjx) 

Cl (a;-xi)(l-a;) 


y e(c) + (1^2) 


£(c) = £fc)l'^68 

_ * * -X2 , a;2S(c) \ , * ( -, , X2 , X2 

— ^162 ■ £3^68 ( 7-7—1“ 7-7 ) + ^i£3 ■ £2^68 I ~1 ~ S(c) + 


+£ 2^3 • £ 1^68 


X — Xi 1 — X 
-X2 X2S(^c) 


1 — X X — Xi 


1 — X X — XiJ ’ 
with S(c) = (xi — X + xi)(l — X + l)/x^ and 

B(r:\ X2 MIq+M^ 


bd 


’(c) 

kt Mt, + Ml-Mt2 


ifesfbd/fl2 — ^) ■ 


d. Vertex functions 7 (rf) and 7 (y^) 


(C19) 


(C20) 


(C21) 


7(d) + 7(/) 


where 


and 



4 ^ 123 ^/ 1 ^ + (1 0 2 ) , 


,-L2 


X2fct 



1 

■ ^t(d) Bt(f) 

1 

to 

_X2k^^ X2k^M'^ 


X2M'^ — M\2 
M^ + M\2-Mt, 





12 



f fbd/fl2 — 1 \ 


(C22) 


^«(d _ -^12 ( f 2 _ 1 \ 
k+ A44 U 


(C23) 


e. Vertex functions 7 (g) ond 7 (i) 


7(3) +7(*) = SlVc^my 


1 + 


(1-x)^ 








't(i) 




fc+A42 


+ 2Y,23^ilj + (lo2) , 

fCo 


(C24) 


where rs^{x) is given in Eqs. (C12|, the ultra-violet regulator is (C13l, rA^(K'*-) = the mass counterterm 

contribution stems from 2g‘^fl‘g = 167 r^)x^, and 


and 



Ml2+M^ 

2M^M‘f2 


(2M^ 


■^ 12 ) 


7^-1 

2A44 


fcaf / fl2 — 1 
2A42(_/V{2 __y)^2^)J ’ 


(C25) 


^«(d _ -^12 / f 2 _ 1 \ 

^3+ W 


(C26) 


Graphs (e) and (/i) result from products of terms S 31 or 7 f 22 with ¥ 31 . The former are independent of the transverse 
momentum, and the latter is odd in transverse momentum. This leads to zero in the integration over k. 
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Appendix D: Three-gluon vertex counterterm 


The ultraviolet regularization dependence in 7 (a) turns out to come only from (see below) 




X2 




(Dl) 


which matches the result found in Ref. |7]. So, 


7(a)div = 8 / dxrst{x)'^^ + (1 ^ 2) , 


1 — X 
X2 


where 


and 


with 


jJciC 2 C 3 [/Ufcj ^ = tT In [Ci2yi2 + CisYis + C23T23] , 


X2 


C 12 = 


Cl3 = 


C 23 — 


>'121 


1 (1 — xY 


1 — X X — Xi X 


X2 


1 (1 — xY 


1 — X X — Xi X 


- 2 


X2 


1 (1 — xY 1 + x"^ 


1 — X 

+ + 


— 

X — Xi X 

xi 

X2 

( 7 : 

to 


S3K12 

1 

5^13 


■ £^1^12 

X2/3 

^23 


■ £Iki2 

1 

Xl/3 


- 2 


The diverging parts are 7(h)div = 7(c)div = 0 and 7(d)div = 27(s)div, as in Ref. [7], with 


(D2) 

(D3) 

(D4) 

(D5) 

(D6) 

(D7) 

(D8) 

(D9) 


7(ff)div =--/VcYi23^5(a;i) / dxrs^[x) 

Jo Jfj.^ 


Tvdhr 






(1-CC)2_ 


+ (1 O 2) . (DIO) 


These results are the same as in Ref. [7] due to the fact that the difference between old and new generators resides 
solely in the RGPEP factors Bt- The new generator leads to BtS that differ from old BtS by the additional terms 
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in numerators, denominators and arguments of exponentials that depend on do not depend on . These 

additional terms do not affect the behavior of BtS when —>■ oo. Terms {i) and (/) are not divergent. 

In summary, in the UV-limit of k'^ —>■ oo, the integrands in all vertex functions behave in the same way as the 
corresponding ones obtained using the old generator [7] . Hence, the divergent part of the vertex counterterm, denoted 
by 7oodiv) is also the same, c.f. Appendix C in Ref. [?]• 

The divergent part of the vertex counterterm is defined by the condition 


7(a)div T i^ 7 ( 5 )div T Toodiv T (1 ^ 2) — 0 . (Dll) 

One thus finds the Hamiltonian vertex counterterm whose vertex function is 

— TT /\ 

7oo = ^123-5- In — [11 + ^(2;i)] + 7finite > (D12) 

3 /i 

where /r denotes the arbitrary separation point between the range of integration over large that extends up to A 
and the finite range of integration where no dependence on A may arise. The finite part of the counterterm, yfinitej 
removes the artificial dependence on /r. The function h{xi) is 


h{xi) = 6 / dxrst{x) 

J X\ 


2 1 1 

1-'-'— 

1 — X X — Xi X 


9 / dxrsfj,{x) 
Jo 


1 

1 — X 


+ (1 O 2) 


(D13) 


The vertex function in the complete vertex counterterm is 


l{j) — 7oo + (1 -t 2) . 


(D14) 


Appendix E: Thrid-order contributions to gt 


The Hamiltonian coupling constant gt is extracted from the term Yt in Eq. (40) that is linear in ^12 i n the limit 
0. Not every term shown in Fig. [^contributes to the running coupling defined this way. Appendix C3 shows 
that terms (e) and (h) do not contribute. Furthermore, terms (6), (c), (/) and (i) vanish faster than linearly in the 
limit k ,-^2 0. Thus, only the terms (a), {g) and (d) contribute to gt- Contribution to gt of each and every term is 





TPSTPPT' 



+ 

8 X (a) + 

2x (b) + 

2 X (c) + 

4 x (d) 


1 X (e) + 

2 x(f) + 

2 x(g) + 

1/2 X (h) 


1 X (i) + 

1 X (j) + 

o(9o) 



FIG. 8: Effective three-gluon vertex (expansion up to third order). 


extracted in three steps: 1) calculation of the coefficient of Kj ^2 ll^® integrand in the limit Kj ~2 —>■ 0; 2) integration 
over K-^; 3) integration over x. Contributions of the counterterm are defined by the subtraction at t = to that is 
described in Secs. HV^ and lrVCl 


1. Contribution to the running coupling from 

The expansion in small k \2 in term (a) concerns the factor 

^t(a) ~ ^to(a) ■■ 1 , 




’■68'^16'^ 


which leads to 


^t(a) ^to(o) I 

kr ' 






(1 — x) Xi 


K 


X2 


12 


Xi 


K{a) ^to{a) 

kr 






(El) 


(E2) 
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This is integrated over yielding 


lim 


Btia) - B, 


to(a) 


k+^ 


^68 


Al(x, Xi) TT f 6 ^^k (2 + - -- In 




where Ai{x,xi) = x'^il — and A 2 {x,xi) = 2x^{l — x)^ Here we have made use of the formula in 

Appendix |E 5[ The running coupling contribution is extracted from 


tio-fO lOTT'^ 2 


j1_ 1 

IGtt^ 2 


fClC2C3 

2 


lim Ct,a(xo,Ki 2 ) , 


(E4) 


where 


TT , t 


lim Ct,a(xi,Ki 2 ) = 7 In — / dxrst(,x) 


(l-x) 


>-0 


4 ^0 Jxi 

Xi{l-X) 

X2X 2^(“) 


X2 


1 - 


Xi{l x)l\jji j lxi(l x) J Jjj 
^ I “r ^ ^12^(a) 


+ (10 2)+ o{Ki2) ■ 


X2X 2 J 2 XX2 

J. 


The contraction of indices of the tensor structure simplifies the limit of Ct^ai^i, k;(‘ 2 ) to 


Ct,a(a;i,Kl^2) ^ jln;^ 
4 to 


f dx rsY ft 2 ix,xi,e-^) + (1 O 2) 

V a^i 


where 


/ I / I \ * * _L * *_L / * *_L / 

i2(^)^i+ ) — Ci2 eie2 £3 ~ ni3 £^£3 ^2 /a;2 — C23 £2£3 £1 /^ci ) 


(E5) 


(E6) 


(E7) 


with the coefficients Cij given in Eqs. (D4|-(D6). Finally, the integration over x leads to the running coupling 
contribution of term (a), 


9t,(a) 


1 t 

9to,{a) + ^§7+2^^ [-ll + 3Xa{a:o)] In— 


(E8) 


with 


Xa{xo) = / dx r^y [2/(1 - a;) + l/(a: - a:o) + 1/a;] + (xo O 1 - a;o) . 


(E9) 


2. Contribution to the running coupling from 'y(^d) 


In this case A^|s “ -^L = M^/x 2 + J^i 2 ' The limit —)■ 0 in Eq. (C22) concerns the factor, 


lim 


(d) ^tp {d) 


rp2 U+2 

*^2^3 


1 - 


^ (l\ ^ E^fl x)"^ f /o/fcd.O ffbd\ 

tmI "Ho? ?ri+i7~J ■ ' ’ 


Integration over k'^ yields: 


d^K'^ lim 


Bt (d) Btg (d) 


0.2^3 


= — a:^(l — a:)^^ In — . 

^ ^ 4 to 


(Ell) 


The last step is the integration over x of 7 (d) — 70 , (d) in this limit. The corresponding contribution to the running 
coupling is 


9t, (d) 


9to,{d) 



(E12) 
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3. Contribution to the running coupling from 7 (g) 


The calculation is analogous to the previous cases. The limit Kj ~2 —>■ 0 concerns the difference of renormalization 
group factors and produces 




>0 kt^ 


2M'^ 


Integration over k'^ of the first two terms gives zero, and the only contributing part is 


lim 


(g) ^tp (g) 

k+^ 


TT , t 


= -x^l-xY-ln 


4 to 


The resulting contribution to the running coupling of term (g) is 


5^ A. 1 , ^ f 11 r , / N [1 111 

9 ,m = «..(.)2‘%r¥+y„ U + 


(E13) 


(E14) 


(E15) 


4. Sum of contributions in App. |E 1[ |E 2| and |E 3| 

Denoting gt^ by go, the sum of contributions (a), (d) and (g) up to order g^ gives 

9t =90- [11 + h{xo)] In ^ . 


(E16) 


5. Useful formula 

Integrals of differences of exponentials that appear in the RGPEP for massless quanta, are evaluated taking advan¬ 
tage of the formula 


= (E17) 

This formula is a consequence of the RGPEP design that secures absence of large perturbative contributions in the 
matrix elements near diagonal of the effective Hamiltonian matrix evaluated in the basis of the Fock space built 
using creation operators for effective particles [2T]. Namely, the arguments of form factors / vanish quadratically as 
functions of the corresponding perturbative denominators. 
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